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A fully relativistic model for the description of exclusive (e, e'p) reactions off nuclear targets at 
high energies and momentum transfers is outlined. It is based on the eikonal approximation for the 
ejectile scattering wave function and a relativistic mean-field approximation to the Walecka model. 
Results for 12 C(e,e'p) and 16 0(e,e'p) differential cross sections and separated structure functions 
are presented for four-momenta in the range 0.8 < Q 2 < 20 (GeV/c) 2 . The regions of applicability 
of the eikonal approximation are studied and observed to be confined to proton knockout in a 
relatively small cone about the momentum transfer. A simple criterium defining the boundaries of 
this cone is determined. The Q 2 dependence of the effect of off-shell ambiguities on the different 
(e, e'p) structure functions is addressed. At sufficiently high values of Q 2 their impact on the cross 
■ sections is illustrated to become practically negligible. It is pointed out that for the whole range of 

Q 2 values studied here, the bulk of the relativistic effects arising from the coupling between the lower 
components in the wave functions, is manifesting itself in the longitudinal-transverse interference 
term. 
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I. INTRODUCTION 
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Exclusive A(e, e'p)B reactions from nuclei constitute an invaluable tool to probe a wide variety of nuclear phenomena. 
At low values of the virtual photon's four-momentum transfer Q 2 = q 2 — lu 2 and, accordingly, larger distance scales, 
the quasi-elastic A(e, e'p) reaction probes the mean-field structure of nuclei. From systematic investigations for a 
large number of target nuclei a richness of precise information about the independent-particle wave functions and 
spectroscopic strengths was assembled jjj. At high Q 2 and decreasing distance scales, the scope of exclusive (e, e'p) 
measurements shifts towards studies of (possible) medium dependencies of the nucleonic properties, and, effects like 
"^T* ! color transparency and the short-range structure of nuclei. Within the context of exclusive (e, e'p) reactions, "color 
transparency" stands for the suggestion that at sufficiently high values of Q 2 the struck proton may interact in an 
^ | anomalously weak manner with the "spectator" nucleons in the target nucleus [0] . 

The extraction of physical information from measured A(e,e'p)B cross sections usually involves some theoretical 
modelling of which the major ingredients are the initial (bound) and final (scattering) proton wave functions and 
the electromagnetic electron- nucleus coupling. At lower values of Q 2 , most theoretical work on (e,e'p) reactions was 
performed in the so-called distorted- wave impulse approximation (DWIA). The idea behind the DWIA approach is 
that the inital (bound) and final (scattering) state of the struck nucleon can be computed in a potential model, 
whereas for the electron-nucleus coupling an "off-shell corrected" electron-proton form can be used. The wealth 
of high-quality (e, e'p) data that electron-scattering experiments have provided over the last 20 years, made sure 
that the DWIA models are well tested against experimental data. For higher values of the energy and momentum 
transfer (Q 2 > 1 (GeV/c) 2 ), most theoretical (e, e'p) work starts from the non-relativistic Glauber theory [||. This 
theory is highly successful in describing small angle proton-nucleus scattering at higher energies [Q and is conceived 
as a baseline for calculating the effect of final-state interactions in high-energy (e, e'p) reactions. Glauber theory 
is a multiple-scattering extension of the standard eikonal approximation that relates through a profile function the 
ejectile's distorted wave function to the elastic proton scattering wave function HH |[. The Glauber method has 
frequently been shown to be reliable in describing A(p,p') processes. Several non-relativistic studies P~|Tl[ have 
formally investigated the applicability of the Glauber model for describing A(e, e'p) reactions at higher energies and 
momentum transfers. Recently, the first high-quality data for exclusive 16 0(e,e'p) cross sections at higher four- 
momentum transfer (Q 2 > l(GeV/c) 2 became available |l2]] . Below, we will compare results of relativistic eikonal 
calculations with these data. We believe that this comparison between model calculations and data provides a stringent 
test of the applicability of the eikonal approximation in describing (e,e'p) reactions.. 

Since relativistic effects are expected to become critical in the GeV energy domain, we explore the possibility of 
developing a fully relativistic model for A(e, e'p) processes, thereby using the eikonal limit to solve the equations for 
the final-state wave functions. We employ a relativistic mean- field approximation to the Walecka model to determine 
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the bound state wave functions and binding energies, as well as nucleon and meson potentials. The same mean-field 
potentials are then also used to compute the scattering wave function in the Dirac eikonal limit. The work presented 
here is a small initial step towards the formulation of a fully microscopic relativistic model for the description of (e, e'p) 
reactions that could possibly bridge the gap between the low and intermediate-energy regime. The model developed 
in this work can be formally applied in a wide Q 2 range. As a matter of fact, we employ the relativistic eikonal 
method to estimate the sensitivity of (e, e'p) observables in the few GcV regime to a number of physical effects, 
including off-shell ambiguities and relativity. We adopt different prescriptions for the electron-nucleus coupling in 
our calculations. By doing this, we estimate the sensitivity of the observables to the theoretical uncertainties that 
surround the choice of the off-shell electron-proton vertex. It is often claimed that off-shell ambiguities decrease in 
importance as the four-momentum transfer increases. Here, we make an attempt to quantify the relative importance 
of the off-shell effects for the (e, e'p) structure functions by comparing results obtained with different off-shell electron- 
proton couplings. Hereby we are primarily concerned with the question how big the uncertainties remain when higher 
and higher four-momentum transfers are probed. 

In Sect. H we introduce a relativistic eikonal formalism for ca lcula ting A(e, e'p) observab les. This includes a 
discussion of the method employed to determine the bound (Sect. [IB ) and scattering (Sect. 1IC) states. Various 
forms for the photon-nucleus interaction vertex are introduced in Sect. |[ID] , where special attention is paid to the issue 
of cu rrent conservation. In Sect. HI we present the results of our 12 C(e, e'p) and 16 0(e, e' p) nu merical calculations. In 
Sect. IIIB we focus on the issue of the Q 2 evolution of the off-shell ambiguities. In Sect. |III C we compare the results 
of a fully relativistic calculation with a calculation in which the explicit coupling between the lower components in 
the inital and final state are neglected. Finally, our concluding remarks are summarized in Sect. [Tv| . 



II. FORMALISM 



A. Reaction observables and kinematics. 



In this work we follow the conventions for the (e, e'p) kinematics and observables introduced by Donnelly and Raskin 
in Ref. jlj|. The four-momenta of the incident and scattered electrons are labeled as K^(e,k) and K ^(e',k'). The 
electron momenta k and k' define the scattering plane. The four-momentum transfer is given by = — K M = 
P%_1 + Pf ~ P^i where P% and P^_ 1 are the four-momenta of the target and residual nucleus, while Pj is the four- 
momentum of the ejected nucleon. Also, g M = (cj, q), where the three-momentum transfer q~k — k' = fcyi-i +kf — kA 
and the energy transfer oj = e — e = Ea-i + Ef — Ea are defined in the standard manner. The xyz coordinate system 
is chosen such that the z-axis lies along the momentum transfer q, the y-axis lies along k x k' and the x-axis lies in 
the scattering plane; the reaction plane is then defined by kf and q. The Bjorken-Drell convention jlj] for the gamma 
matrices and Dirac spinors is followed, so that the normalization condition for Dirac plane waves, characterized by a 
four-momentum and spin-state , is u(K>>, , S^) = 1. 

In the one-photon-exchange approximation, the process in which a longitudinally polarized electron with hclicity h, 
impinges on a nucleus and induces the knockout of a single nucleon, leaving the residual nucleus in a certain discrete 
state, can be written in the following form Ea] : 
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with Of the angle between kf and q, and gm the Mott cross section 
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with 9 e the angle between the incident and the scattered electron. The electron kinematics is contained in the 
kincmatical factors 
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whereas the structure functions are defined in a standard fashion as 

Kl = \ P (q)M 2 , (10) 

7e T = |J(g;+l) /i | 2 + |J(g;-l) /i | 2 , (11) 

TItt =2Re{J*{q;+l) fi J(q;-l) fi } , (12) 

Ktl = -2Re{p*(q) fi (J(q;+l)fi - J® -1))/;} , (13) 

ftr, = |J(<f; +1)/;] 2 - I J(g; -1)/*| 2 , (14) 

^tl' = -2Re{p*(q) fi {J{q;+l) fi + J(g;-1))/;} , (15) 

where p(q)fi is the transition charge density, while J(q;m — ±l)/j is the transition three-current expanded in terms 
of the standard spherical components. 



B. Bound state wave functions. 



A relativistic quantum field theory for nucleons (ip) interacting with scalar mesons (<fi) through a Yukawa coupling 
t/jipcj) and with neutral vector mesons (Vu) that couple to the conserved baryon current ip^f^ip, can be described through 
a lagrangian density of the type fl^Juj 

£ = $(t/B - M)i> + i(a„ w - m 2 j 2 ) - \g^g» v 

+ \™lV»V» - frftyrfV + g s ^U , (16) 

with M, m s and m v the nucleon, scalar meson and vector meson masses, respectively, and = d^V" — d v V tl the 
vector meson field strength. The scalar and vector fields may be associated with the a and lo mesons. The model can 
be extended to include also tt and p mesons, as well as the coupling to the photon field. The corresponding lagrangian 
has the form 

C = Co + ~{d„w ■ d^TT - mln ■ tt) - ig^ l5 r ■ 7n/> - -B^ ■ B^ 
+ \m 2 p b p ■ - ~g p jjw ■ - ~ V" 

-eA^j^(l + r 3 )V + (K x B^) 3 + (tt x (0"tt + 5p (vr x 6"))) 3 ] • (17) 

Here tt, b^, A p , F M „ are the pion, rho, Maxwell and electromagnetic fields. Further, B^ v = d^b v — d v b^ — g p (b^ x b v ) 
is the p-meson field. 

At sufficiently high densities, the meson field operators can be approximated by their expectation values. Within 
the context of the relativistic Hartree approximation, it can be shown that when starting from the langrangian ( p"^ ) 
the following Dirac equation for the baryon field \& results [jl6| : 
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fa^ — M — Yj H ] * = 



where the self-energy is defined as 



(18) 
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Assuming that the nuclear ground state is spherically symmetric and a parity cigcnstate, it can be shown that the 
pion field does not enter in the Hartree approximation. Furthermore, the meson fields only depend on the radius, and 
only the time component of the vector fields contribute. The time-independent Dirac equation can then be written 
as : 



-ia-V + g v V°(r) + ±g p T a b 0a (r) 
+ ±e(l + T 3 )A (r) +1 (M~g s ^(r)) 



The general solutions to a Dirac equation with spherically symmetric potentials have the form 
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where n denotes the principal, n and m the generalized angular momentum and t the isospin quantum numbers. The 
y± Km are the well-known spin spherical harmonics and determine the angular and spin parts of the wavefunction, 

y Km = ^2 < lmi^m s \lpm > Y limi xi ms , 
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The Hartree approximation leads to the following set of coupled equations for the different fields m& 
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The above equations constitute the basis of the relativistic mean- field approach to the lagrangian of Eq. (|17|) . 

A new computer program to solve the above set of coupled non-linear differential equations was developed. Starting 
from an initial guess of the Woods-Saxon form for the scalar and vector potential, the Dirac equations can be solved 
iteratively using a shooting point method. Analytic solutions to the equations in the regions of large and small r allow 
to impose the proper boundary conditions. Once the nucleon wave functions are obtained, the densities and meson 
fields can be re-evaluated. This procedure is repeated a number of times until convergence for the energy eigenvalues 
is reached. We adopt the values for the a, u> and p masses and coupling constants as they were introduced by Horowitz 
and Serot |j§. 

For the ^} and 16 nuclei, the newly developed C-code SDR performed all integrations for a radial extension of 
the nucleus of 20 fm and a stepsize of 0.01 fm. The coupled Dirac equations were solved for a shooting point lying at 
2 fm using a fourth order Runge-Kutta algorithm. As a convergence criterium we imposed a tolerance level as small 
as 0.001 MeV on all single-particle energy levels. The computed densities for the nuclei 12 C and 16 0, are depicted in 
Fig. |l|. We have verified that these results are comparable to those produced by the TIMDRA code [^6|, which is widely 
used to solve the set of Eqs. (p^). 

C. The eikonal final state. 

To construct the scattering states for the ejected nucleons, we consider the hamiltonian ( |20| ) that was already used 
to calculate the bound state wave functions 

H = -ia- V + 7°M + 7°£#(r) , (24) 

where the self-energy Y,h(t) is given by 

E H (r) = -g s Mr)+9vloV°(r) + l -g plo r a b 0a (r) + ±e 7 o(l + r 3 )A°(r) . (25) 
With the formal substitutions 

v s(r) = -9s<Po , 

V v (r) = g v V Q (r) + I 5p 6 ( r )(-1)*«-V2 + e A (r)(t a + ±) , (26) 



the time independent Dirac equation for a projectile with relativistic energy E = Vfc 2 + M 2 and spin state s, can be 
cast in the form 

=[a-p + 0M + pV s (r) + V v {r)W ] , (27) 
where we have introduced the notation 6^ for the unbound Dirac states. The computed scalar and vector potentials 

k,s 

for the 12 C and 16 O nuclei are displayed in Fig. || 

After some straightforward manipulations, a Schrodinger-like equation for the upper component can be obtained 



7^+V c + V so (a-L-if-p) 



u (+) = *l u (+) , (28) 

k,s 2M k,s ' y ' 



where the central and spin orbit potentials V c and V so are defined as 

V c (r)=V s (r) + *V v (r) + V ^ 2 - V ^ 2 , 
M 2M 

2M [E + M + V s (r) — V v (r)J r dr 

In computing the scattering wave functions, we use the scalar and vector potentials as obtained from the iterative 
bound state calculations. As a result the initial and final state wave functions are orthogonalized and no spurious 
contributions can be expected to enter the calculated cross sections. 
Since the lower component is related to the upper one through 

wi +) = 8 ■ pui +) , (30) 

*.« E + M + V.-Vv*'' 



5 



the solutions to the equation (E8h determine the complete relativistic ei gen value problem. So far no approximations 
have been made. Various groups |jl7|-|l9| have solved the Dirac equation ( |28|) fo r the final scattering state using Dirac 
optical potentials derived from global fits to elastic proton scattering data ]2C| . Not only are global parametrizations 
of Dirac optical potentials usually restricted to proton kinetic energies T p < 1 GeV, calculations based on exact 
solutions of the Dirac equation frequently become impractical at higher energies. This is particularly the case for 
approaches that rely on partial- wave expansions in determining the transition matrix elements. To overcome these 
complications, we solve the Dirac equation ( p8f ) in the eikonal limit pl[p2| . In intermediate-energy proton scattering 
(T p « 500 MeV) the eikonal approximation was shown to reproduce fairly well the exact Dirac partial wave results 
p3| . Following the discussion of Ref. p3| , we define the average momentum K and the momentum transfer q in terms 
of the projected initial (ki) and final momentum (kf) of the ejectile 



q = k f - ki , 
K=\(k S +q) 



(31) 
(32) 



In the eikonal, or, equivalently, the small-angle approximation (g> ki) the following operatorial substitution is made 
in computing the scattering wave function 

p 2 = [(p-K) + K} 2 — >2K-p-K 2 . (33) 
After introducing this approximate relation, the Dirac equation for the upper component (|28| ) becomes 

[-lit ■ V - K 2 + M(V C + V so [a -(rxit)-if- = , (34) 



where the momentum operators in the spin orbit and Darwin terms are substituted by K. Remark that the above 
equation is now linear in the momentum operator. In the eikonal limit, the scattering wave functions take on the 
form 



(35) 



Inserting this into Eq. (|34|), yields an expression for the eikonal phase [|TJ. Defining the z-axis along the direction of 
the average momentum K, this phase can be written in an integral form as : 



tS(b,z) = -%^ I ,):'T,{h. :') \ :.,U>. :')[*■ [bx K ) ■ ■ ,AV | . 



(36) 



where we have introduced the notation r= (b, z). The scattering wave function, which is proportional to 



k.s 



is normalized such that 
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This wave function differs from the plane-wave solution in two respects. First, the lower component exhibits the 
dynamical enhancement due to the combination of the scalar and vector potentials. Second, the eikonal phase e lS ^ 
accounts for the interactions that the struck nucleon undergoes in its way out of the nucleus. The calculation of the 
eikonal phase (|36| ) involves a transformation to a reference frame other than the usual laboratory or center-of-mass 
frame, namely the frame where the average momentum is pointing along the z-axis. As the eikonal phase has to be re- 
evaluated for every (&, z) point in space, the Dirac eikonal (e, e'p) calculations are very demanding as far as computing 
power is concerned. In evaluating the matrix elements, the radial integrations were performed on a 0.1 fm mesh. It is 
worth remarking that the standard Glauber approach followed in many studies involves an extra approximation apart 
from the ones discussed above. Indeed, in evaluating the eikonal phase from Eq. ( |36"| ) one frequently approximates 
the z-dependence of the potentials by a delta function. 
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D. Off-shell electron-proton coupling 



We express the matrix elements of the nucleon current in the usual form 

< P f S f \J»\PiSi >= UfT^P^P^u, , (39) 

where T M is the electromagnetic vertex function for the nucleon and itj (it/) the nucleon spinors. As discussed in 
many works p^-|2^|, some arbitrariness, often referred to as the "off-shell ambiguity", surrounds the choice for the 
functional form of the vertex function r M . For a free nucleon, can be expressed in several fully equivalent forms 

T^ cl = G M (Q 2 )r - ^F 2 (Q 2 )(P? + Pf) , (40) 

Kc2 = *i(QV + t JL F F 2 {.Q 2 )a^q p , (41) 



2M 

C 3 = ^i(Q 2 )(Pf + + (42) 



where F\ is the Dirac, F 2 the Pauli form factor and k is the anomalous magnetic moment. The relation with the Sachs 
electric and magnetic form factors is established through Ge = Fx — tkF 2 and Gm = Fx + kF 2 , with r = Q 2 /Am 2 . 

When considering bound (or, "off-shell" ) nucleons, however, the above vertex functions can no longer be guaranteed 
to produce the same results. As a matter of fact, explicit current conservation is rather an exception than a rule in 
most calculations that deal with (e, e'p) reactions from finite nuclei. In nuclear physics, the most widely used procedure 
to "effectively" restore current conservation is based on modifying the longitudinal component of the nuclear vector 
current using the substitution 

J z -» -Jo • (43) 

q 

This procedure is partly inspired on the observation that meson-exchange and isobar terms enter the charge current 
operator in a higher relativistic order than they used to do for the vector current. There exist several other prescriptions 
which are meant to restore current conservation. Along similar lines, the charge operator can be replaced by 

Jo -» -J, ■ (44) 

One can also construct a vertex function that garantuees current conservation for any initial and final nucleon state. 
This can be achieved for example by adding an extra term to the vertex P9] 



Kon = FxiQ 2 )^ +i ^F 2 (Q 2 )a^q v + Fl {Q 2 )^ , (45, 



which is also equivalent to the Eqs. ( |40| - [42| ) in the free nucleon case. An operator derived from the generalized 
Ward-Takahashi identity reads pq| 



r^r = V - *2^(Q 2 K^ + [Fi(Q 2 ) i] W + 2 Q V . ,4<i. 



III. RESULTS 



A. Final state interactions and the eikonal approximation 

We start our (e, e'p) investigations within the relativistic eikonal approximation for the kinematics of an 16 0(e, e'p) 
experiment that was recently performed at Jefferson Lab Jl2| . In this experiment, the separated 16 0(e, e'p) structure 
functions are measured at Q 2 — 0.8 (GeV/c) 2 and u> = 0.439 GeV for missing (or, initial) proton momenta p m =\ kf—q \ 
below 355 MeV/c. The variation in missing momentum was achieved by varying the detection angle of the ejected 
proton with respect to the direction of the momentum transfer ("quasi-perpendicular kinematics"). The measured 
cross sections for knockout from the lpx/2 and levels are depicted in Fig. [| along with the predictions of our 

relativistic eikonal calculations. A spectroscopic factor of 0.6 was adopted for all bound levels, and the standard 
dipole form was used for the electromagnetic form factors. At low missing momenta, the eikonal results shown in Fig. 
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|3| produce a fair description of the data. As a comparison, the results of a relativistic plane wave calculation in the 
impulse approximation (RPWIA) are also displayed. Through comparing the plane- wave and the eikonal calculations, 
thereby keeping all other ingredients of the calculations identical, one can evaluate how the eikonal method deals with 
final state interactions (FSI). In the eikonal calculations, the dips of the RPWIA calculations are filled in, and, at low 
missing momenta the RPWIA cross sections are reduced. These two features reflect nothing but the usual impact 
of the final-state interactions on the A(e, e'p) angular cross sections. The limitations of the eikonal approximation 
(q ^> ki) are immediately visible at higher missing momenta (p m > 250 MeV/c). Here, the eikonal cross sections 
largely overshoot both the RPWIA results and the data and should by no means be considered as realistic. It is worth 
remarking that the data closely follow the trend set by the RPWIA curves. As a matter of fact, whereas the eikonal 
calculations predict huge effects from final-state interactions at large transverse missing momenta, the data seem to 
suggest rather the opposite effect. We consider this observation as one of the major findings of this work. 

One may wonder whether the observed behaviour of the eikonal results at higher missing momenta in Fig. || is a 
mere consequence of the small- angle approximation contained in Eq. (|33|), or whether the adopted model assumptions 
for computing the scattering states is also (partly) at the origin of this pathological behavior. To address this question, 
we have performed calculations for various fixed recoil angles 9 defined as 

n Pm ' Q / ah\ 

cosfl = (47) 

\Pm\ m 

The results are displayed in terms of the reduced cross section p which is defined in the standard fashion as the 
differential cross section, divided by a kinematical factor times the "CC1" off-shell electron-nucleon cross section 
of Ref. Q . For the results of Figure ||| we considered in-plane kinematics at a fixed value of the outgoing proton 
momentum (fc/=l (GeV)/c) and an initial electron energy of 2.4 GeV. The variation in missing momentum is achieved 
by changing the q. For recoil angles 9 — 0° ("parallel kinematics") the eikonal calculations do not exhibit an unrealistic 
behavior up to p m =0.5 GeV/c, which is the highest missing momentum considered here. With increasing recoil 
angles, and consequently, growing "transverse" components in the missing momenta the "unrealistic" behaviour of 
the eikonal results becomes manifest. Accordingly, the accuracy of the eikonal method based on the small-angle 
approximation of Eq. ( p3| ) can only be guaranteed for proton knockout in a small cone about the momentum transfer. 
A similar quantitative behaviour as a function of the recoil angle to what is observed in Fig. ^ was reported in Ref. 
H for d(e, e'p)n cross sections determined in a Glauber framework. We conclude this section with remarking that 
the eikonal method does not exclude situations with high initial (or, missing) momenta, it only requires that the 
perpendicular component of ejectiles's momentum kf is sufficiently small. It speaks for itself that such conditions are 
best fulfilled as one approaches parallel kinematics. This observation puts serious constraints on the applicability of the 
Glauber method, that is based on the eikonal approximation, for modelling the final-state interactions in high-energy 
(e, e'p) reactions from nuclei. However, it should be noted that our framework does use purely real scalar and vector 
potentials. More realistic scattering potentials demand an imaginary part that accounts for the inelastic channels that 
are open during the reaction process. The Glauber approach effectively includes these inelastic channels and on these 
grounds one may expect that its range of applicability is somewhat wider than what is observed here. With the eye on 
defining the region of validity for the eikonal approximation more clearly, we have studied differential cross sections 
for various Q 2 . In Fig. [|, we display the computed differential cross sections for the 12 C(e,e'p) 11 B(lp~^ 2 ) process 

against the missing momentum for Q 2 varying between I and 20 (GeV/c) 2 . Hereby, quasi-elastic conditions were 
imposed. The arrow indicates the missing momentum where the slope of the eikonal differential cross section starts 
deviating from the trend set by the RPWIA cross section. In the light of the conclusions drawn from the comparison 
between data and the eikonal curves in Fig. ^, the eikonal results should be regarded with care beyond this missing 
momentum. Furthermore, it is clear that the change in the slope of the angular cross section becomes more and more 
pronounced as Q 2 increases. It is apparent from Fig. |5| that the eikonal differential cross section changes slope at 
about p m =250 MeV/c for all values of Q 2 considered. We remark that we imposed quasi-elastic conditions for all 
cases contained in Fig. |. As a consequence, the momentum of the ejected nucleon varies quite dramatically as one 
moves up in Q 2 . The uniform behaviour of all curves contained in Fig. fallows one to write down a relation between 
the transferred momentum q and the polar scattering angle 9 : \q\ 9 < 250 MeV rad. This simple relation could serve 
as a conservative guideline to determine the opening angle of the cone in which the outgoing proton momentum has to 
reside to ascertain that the eikonal approximation produces "realistic" results. This limitation of the eikonal method 
can also be inferred from the results contained in Refs. PJB1|]. Indeed, in Figs. 3 and 4 of Ref. Q one can confirm that 
the above relation between \q\ and 9 defines the missing momentum at which a sudden change in the p m dependence 
of the calculated cross sections is observed. The above relation can be understood as follows. In quasi-perpendicular 
and quasi-elastic kinematics, the missing momentum roughly equals the transverse momentum of the ejected nucleon. 
With increasing momentum transfer, the longitudinal momentum of the escaping nucleon increases correspondingly 
while its transverse momentum has to stay smaller than the suggested value of 250 MeV/c. Hence, the sine of the 
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angle between the transferred momentum and the ejectile's momentum has to decrease. Since we are dealing with 
small angles, sin(0) can be approximated by 9. The opening angle of the cone in which the eikonal approximation 
is valid, can be inferred to be independent of Q 2 in the Lorentz frame where the ejected nucleon is at rest. When 
transforming back to the labframe, lateral dimensions become dilated, and, thus, angles contracted. 

B. The Q 2 evolution of off-shell effects 

A major point of concern in any A(e, e'p)B calculation are the ambiguities regarding the off-shell electron-proton 
coupling. Most calculations do not obey current conservation and a variety of prescriptions have been proposed to 
partially cure this deficiency. Here we adopt a heuristic view and estimate the sensitivity of the calculated observables 
by comparing the results obtained with different viable prescriptions for the electron-proton coupling. Amongst 
the infinite number of possible prescriptions for the off-shell electron-proton coupling we have selected four that are 
frequently used in literature. Figure ^ shows the separated structure functions for lpin knockout in the kinematics 
of Fig. [}]. Current conservation was imposed by either modifying the longitudinal component of the vector current 
operator (hereafter denoted as the "JO method"), or by modifying the charge operator (hereafter denoted as the "J3 
method"), along the lines of Eqs. ( ^3| ) and (fli]). Note that for the operator of Eq. (|45|), both methods yield the same 
results, since, by construction, this operator is current conserving, regardless of the method adopted to compute the 
wave function for the initial and final state. 

Turning to the results shown in Fig. ^|, the predicted strengths in the longitudinal structure functions Rl and 
Rtl depend heavily on the choice made for the electron-proton coupling. For the CC1 prescription, for example, 
the values obtained with the J3 method are several times bigger than those obtained within the JO method. The 
predicted differences among the various current operators within one scheme ("JO" or "J3") are also sizeable. The 
ambiguities are, however, much smaller for the calculations performed with the J z — > ^ Jo substitution. This clearly 
speaks in favor of this recipe which is mostly used in A(e,e'p) calculations. The Rtt and Rt structure functions are, 
obviously, insensitive to whether the "JO" or "J3" method is adopted. All adopted electron-proton couplings but the 
CCl one produce the same results in the Rt and Rtt responses. 

With increasing Q 2 and the corresponding decreasing distance scale, the role of off-shell ambiguities in the photon- 
nucleus coupling is expected to decline and the impulse approximation is believed to become increasingly accurate. 
In order to investigate the degree and rate to which this virtue may be realized, we have performed calculations for 
kinematics in the range of 0.8 < Q 2 < 20 (GeV/c) 2 . We use two techniques to estimate the relative importance of the 
off-shell effects as a function of Q 2 . First, results computed with the "JO" and "J3" method can be compared. Second, 
predictions with various choices for the electron-proton coupling are confronted with one another. The validity of the 
IA is then established whenever the final result happens to become independent of the adopted choice. In order to 
assess the degree to which this independence is realized, we have considered ratios of structure functions for some fixed 
kinematics but calculated with different choices for the electron-proton coupling. As a benchmark calculation, we 
have computed 12 C(e, e'p) 11 B(lp~^, 1 2 ) observables in quasi-elastic kinematics for several values of the four-momentum 
transfer. The results are shown in Figs. Q and |[ Fig. shows for several observables the ratio of the values obtained 
with the "J3" scheme to the corresponding prediction using the "JO" scheme. Fig. || shows the ratio of the strengths 
obtained with the CCl vertex function compared to the corresponding predictions with the CC2 form. Remark that 
in the limit of vanishing off-shell effects, these ratios should equal one. It is indeed found that the calculations that are 
based on the substitution J z — ► ^ Jo, tend to converge to those based on Jo — > £ J z with increasing energy transfer. 
This is particularly the case at low missing momenta, where the decrease in the longitudinal response is almost 
exponential. The overall behaviour is identical for the higher missing momentum case (p m — 150 MeV/c), but the 
rate of decrease is somewhat slower. This can be attributed to the fact that at higher momenta, hence, greater angles, 
the transverse components of the vertex functions play a more important role. Looking at Fig. ^| one can essentially 
draw the same conclusions. The predictions with the different prescriptions also converge to each other as the energy 
is increased. Again this convergence is more pronounced for the low missing momentum case. This feature is most 
apparent in the purely transverse channel, which dominates the cross section at sufficiently high energies. It appears 
thus as if off-shell ambiguities, speaking in terms of strenghts and absolute cross sections, are of far less concern at 
higher Q 2 than they used to be in the Q 2 < 1 (GeV/c) 2 region, where most of the data have been accumulated up 
to now. The interference structure functions Rtt and Rtl are subject to off-shell ambiguities that are apparently 
extending to the highest four-momentum transfers considered here. This feature was already established in Ref. |29f | 
and explained by referring to the large weight of the negative energy solutions in the interference structure functions 
Rtl and Rtt- 







C. Relativistic Effects 



Recently, there have been several claims for strong indications for genuine (or, "dynamic") relativistic effects 



in A(e,e'p) observables [12 32 -3J]. In an attempt to implement some of these effects in calculations based on a 
Schrodinger picture, several techniques to obtain a "relativized version" of the electron-nucleus vertex have been 
developed. In leading order in a p/M expansion these "relativized" electron-nucleus vertices typically miss the 
coupling between the lower components in the bound and scattering states. For that reason, we interpret the effect 
of the coupling between the lower components in the bound and scattering states as a measure for the importance of 
relativistic effects. In Fig. || we display results of fully relativistic 12 C(e, e'p) 11 B(lp~y 1 2 ) calculations and calculations 
in which the specific coupling between the lower components in the bound and scattering states have been left out. 
We consider quasi-elastic conditions and study the Q 2 evolution of the structure functions for two values of the 
missing momentum (p m = and 150 MeV/c) both corresponding with small recoil angles. Hence, the results of Fig. O 
refer to kinematic conditions for which the eikonal approximation is justified. A rather complex and oscillatory Q 
dependence of the relativistic effects emerges from our numerical calculations. Looking first at the p m ss MeV/c 
case, which nearly corresponds with parallel kinematics, we observe that for both the longitudinal and transverse 
structure functions, the impact of the coupling amongst the lower components first increases, and then tends to 
become fairly constant for higher values of lu. The genuine relativistic effect stemming from the coupling between 
the lower components in the initial and final states is larger in the longitudinal than in the transverse channel. It is 
noteworthy that in the cross section the impact of the "relativistic dynamical effects" never exceeds the 10 % level. If 
we turn our attention to the interference structure functions Rtl and Rtt, the relativistic effects grow in importance. 
Especially for the Rtl structure function the effects are large and extend to the smallest values of Q 2 considered here. 
This enhanced sensitivity of the Rtl response to relativistic effects, even when relatively low values of Q 2 are probed, 
complies with the conclusions drawn in other studies |28| , ^2] , ^5| -[37| . Also the tendency of the relativistic effects to 
increase the cross section when higher values of p rn are probed complies with the findings of earlier studies p§| ] . A 
quantity that is relatively easy to access experimentally and depends heavily upon the Rtl term, is the so-called 
left-right asymmetry Alt 

= <r(0 = 0°) - = 180°) = VTLRTL , , 

LT a(<j) = 0°) + a{(j) = 180°) v L R L + v t Rt + v TT Rtt ' 

In Fig. [l0]we have plotted the left-right asymmetry for both lpi/2 and lp^/2 knockout from 16 in the kinematics of 
Fig. ^. It is indeed verified that the asymmetry is very sensitive to relativistic effects. As has been reported, relativistic 
effects enhance the asymmetry further, and this enhancement is more pronounced for the lpi/2 knockout reaction. 
The role played by the lower components in this dynamical enhancement of the left-right asymmetry can be further 
clearified by looking at the results of Fig. [ll| In this figure, we plot the left-right asymmetry for lp^/2 knockout 
from 12 C, for different Q 2 and quasi-elastic conditions. Looking at the fully relativistic curves, we observe a gradual 
decrease of the asymmetry with increasing Q 2 . At the same time, the relative contribution of the "non-relativistic" 
contribution to Alt diminishes. This feauture indicates that the asymmetry Alt is nearly exclusively generated by 
the coupling between the lower components as Q 2 



IV. SUMMARY 



We have outlined a fully relativistic eikonal framework for calculating cross sections for (e, e'p) reactions from 
spherical nuclei at intermediate and high four-momentum transfers and carried out 12 G(e,e'p) and le O(e,e'p) calcu- 
lations for a variety of kinematical conditions, thereby covering four-momentum transfers in the range 0.8 < Q 2 < 20 
(GeV/c) 2 . Our results illustrate that the validity of the eikonal method is confined to proton emission in a cone with 
a relatively small opening angle about the direction of the virtual photon's momentum. This observation puts serious 
constraints on the applicability of the Glauber method, that is based on the eikonal approximation, for modelling 
the final-state interactions in high-energy (e, e'p) reactions from nuclei. Incorporation of the inelastic channels in the 
eikonal method is however needed to fully appreciate the limits of the Glauber model, and work along these lines 
is in progress. In line with the expectations, our investigations illustrate that the uncertainties induced by off-shell 
ambiguities on the calculated observables diminish as Q 2 increases. Nevertheless, in the relativistic eikonal framework 
four-momentum transfers of the order 5 (GeV/c) 2 appear necessary to assure that the effect of the off-shell ambiguities 
can be brought down to the percent level. Our theoretical framework permits to assess the impact of the relativistic 
effects over a wide energy range. The impact of the lower components on the (e, e'p) observables is observed to be 
significant over the whole Q 2 range studied. Especially the left-right asymmetry lends itself very well to study these 
effects of genuine relativistic origin. 
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FIG. 1. The calculated scalar (p s ), baryon (ps), rho (P3) and proton (pp) density distributions in 12 C and 16 0. 

FIG. 2. The radial dependence of the scalar and vector potentials (in absolute values) as obtained from relativistic Hartree 
calculations for 12 C and ls O. 

FIG. 3. Measured le O(e,e'p) cross sections compared to relativistic eikonal and RPWIA calculations at e = 2.4 GeV, q = 
1 GeV/c and ui = 0.439 GeV in quasi-perpendicular kinematics. The calculations use the current operator CC1. The data are 
from Ref. jjj. 

FIG. 4. The reduced cross section for the 16 0(e, e'p) 15 N(lp~^ 2 ) reaction versus missing momentum at three values of the 
recoil angle 9. A fixed outgoing proton momentum of |p| = 1 GeV was considered. The solid line shows the fully relativistic 
eikonal calculation, while the dashed one shows the RPWIA results. The calculations use the CC1 prescription. 

FIG. 5. The differential cross section for the 12 C(c, e'p) 11 B(lp^, 1 2 ) reaction versus missing momentum at six different values 
for Q 2 . Quasi-elastic conditions and perpendicular kinematics were considered. 

FIG. 6. The different structure functions versus missing momentum for lpi/2 knockout from 16 O in the kinematics of Fig. [| 
The calculations in the left column imposed current conservation by replacing the longitudinal component of the vector current 
operator (Eq. [|3]), while for the results in the right column the charge density operator was modified according to Eq. (^ij). 
The curves refer to the different off-shell prescriptions as they were introduced in Sect. 



II D 



FIG. 7. The Q 2 dependence of the sensitivity of the calculated (e, e'p) structure functions to the choice for the electron-nucleus 
vertex for IP3/2 knockout from 12 C in quasi-elastic kinematics. The curves show for the various observables the ratio of 
the predictions with the "J3" method to those obtained with the "JO" method. Solid (dashed) line corresponds with the 
p m = MeV/c (p m = 150 MeV/c) situation. 

FIG. 8. The Q 2 dependence of the sensitivity of the (e, e'p) structure functions to the choice for the electron-nucleus vertex 
for IP3/2 knockout from 12 C. The curves display the ratio of the predictions using the vertex function r£ cl to those using T[f c2 . 
Solid (dashed) lines correspond with p m = MeV/c (p m = 150 MeV/c). 

FIG. 9. The Q 2 dependence of the sensitivity of the (e, e'p) structure functions to dynamical relativistic effects. The curves 
show for IP3/2 knockout from 12 C the ratio of the fully relativistic results to the predictions when the coupling between the 
lower components is neglected. The solid (dashed) line presents results for the situation p m = MeV/c (p m = 150 MeV/c). 

FIG. 10. The left-right asymmetry Alt for both lpi/2 and IP3/2 knockout from le O in the kinematics of Fig. ^. The data 
are from Ref. [111. 



FIG. 11. The left-right asymmetry Alt for IP3/2 knockout from 12 C for different Q 2 , under quasi-elastic conditions and 
perpendicular kinematics. 
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